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ABSTRACT 

We present results of an imaging observation of the central region of a giant radio 
galaxy B1358+305. The classical, standard scenario of Fanaroff-Riley II radio galaxies 
suggests that shock produced hot electrons contained in a radio galaxy are a good 
reservoir of the jet-supplied energy from active nuclei. The aim of our observation is to 
search for the Sunyaev-Zel'dovich effect induced by these hot electrons. The observation 
was performed at 21 GHz with the Nobeyama 45-m telescope. Deep imaging observation 
of a wide region of size 6.7' x 6.7' with the beam size ^hpbw = 81.2" enables the 
most detailed examination of the possible thermal energy of electrons contained in a 
radio galaxy. The resultant intensity fluctuation is 0.56 mJy beam~^ (in terms of the 
Compton y-parameter, y = 1.04 x 10~^) at a 95 % confidence level. The intensity 
fluctuation obtained with imaging analysis sets the most stringent upper limit on the 
fluctuations in the central region of a giant radio galaxy obtained so far, and our results 
will be a toehold for future plans of SZE observation in a radio galaxy. 
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Introduction 



Extended lobes of radio galaxies are interesting as a probe of the e nergetics of active galax ies. In 
the standard model of radio galaxies of Fanaroff- Riley type II (FR-II: iFanaroff &: Rilevlll974l ). radio 
lobes are forme d by shock interactions of the jets with the surrounding intergalactic medium at the 
jet extremities (|Scheueiill974l : lBlandford &: Reeall974l ). In this scenario, radio lobes consist mainly 
of the jet-supplied matter that passed the termination shock. The gas density in the radio lobes is 
sufficiently low that radiative cooling is ineffective, resulting in the energy supplied to the radio lobes 
to be well conserved for the lifetime of the radio galaxy. The lobe has a higher pressure than the 
ambient intergalactic medium (IGM) due to sho ck-compression, and it can expand supersonically 
into the ambient IGM ( "overpressured cocoon" : iBlandford &: Reeslll974l : iBegelman &: Cioff]||l989l : 
Nathlll995l : iHeinz. Reynolds. &: Begelmanlll998l : lYamada et al.lll999l ). It is expected to form a shell 
of IGM matter compressed by the external shock exterior to the radio lobes. 

Detections of inverse Compton scattered photons in the X-ray energy bands of radio galaxies 
can drastically improve our understanding of the properties of hotspots and radio lobes. The 
distribution of energy among the magnetic field component and in the non-thermal electrons 
has been investigated by e xamining the spectral energy distribution from the radio to X -ray 



bands ( e.g., for Cygnus A. [Harris. Carilli. fc Per 



3C123, Hardcastle. Birkinshaw. and Worral| 



Wilson. Young. &: 



ShoDbeli 



2001 



2001 



1994 : IWilson. Young. &: Shoobelll l2000l : for 



for 3C295, 



Harris et al 



2001 



B; for Pictor A, 



; for 3C120. lHarris et al.lll999l : for 3C390.3. lHarris. Leighlv. k Leahy 
19981 ). The results of these studies revealed that high energy elect rons of 7 ~ 10^ — 10^ have an en- 



ergy com£arable_toor^ field (e.g.. lTashiro et al.lll998l : lGrandi et al 



2003 



Croston et al 



20051; 



Isobe et al. 



2006 



2OO9I ). If the energy spectrum of non-thermal electrons 



is given by a power-law {dN[E)/dE oc E~^), most of the non-thermal electron energy is carried by 
those electrons characterized by energies near the lower cut-off energy of N(E) for rj > 1. However, 
direct measurements of the thermal and/or low energy (7 ~ 10'^ or less) electrons are quite difficult, 
since the plasma density in a radio galaxy is too low to emit detectable radiation. The standard 
overpressured cocoon models of FR-II radio galaxies are characterized by two types of shocks, i.e., 
the external and the internal (jet-terminal) shock, where the latter is believed to be observed as 
hotspots. In this model, radio galaxies emit synchrotron radiation by shock accelerated electrons 
with an acceleration efficiency = f4yn/tAoti which is inferred to be low (t/syn and C/tot are energy 
of synchrotron-emitting electrons and total internal energy, respectively) . If we adopt the standard 
model and assume the low acceleration efficiency in both shocks (low ^g), we can expect that a large 
amount of thermal electron energy, as well as non-thermal, synchrotron-emitting electron energy. 
In order to investigate the energetics of active galaxies, it is quite important to measure all of the 
energy in the electron component, including those that emitting strongly. 

In this paper we employ the Sunyaev-Zel'dovich effect (SZE) as a tool to measure the energy 
of the electrons in a radio galaxy. The SZE represents the spectral deformation of the Cosmic 
Microwave Background (C MB) radiation due to the in verse Compton scattering of these photons 
by the energetic electrons (jZel'dovich &: Sunyaevlll969l ) in the galaxy. The intensity change of the 
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thermal SZE is classically described as follows (jZerdovich &: Sunvaevlll969l ): 



xe 



y 



(e^ - 1) 



+ 1 
- 1 



(1) 

(2) 



where x = hu/ksTr is the non-dimensional frequency, Tr is the temperature of the CMB, A;^ is 
the Boltzmann constant, Tg is the electron kinetic temperature, ax is the cross section of Thomson 
scattering, Ue is the electron number density, rrie is the electron mass, and pe is the electron thermal 
pressure, respectively, integrated along the line of sight. Equation ([2]) shows that the Compton 
parameter y is proportional to the sum of the thermal pressure of electrons along the line of sight. 
Therefore we can estimate the thermal energy of the electrons contained in a radio galaxy with 
the thermal SZE. Similarly, the decrease in CMB intensity by SZE in the Rayleigh-Jeans regime 
reflects all of the electron energy, not only the high energy electrons which generate X-ray photons, 
but also the lower energy electrons and the thermal electrons. Although the energy distribution of 
the electrons is lost in the SZE (Eq.[2]), it measures the total energy deposited in the radio galaxy. 
In this paper, we ignore the non-thermal and kinetic SZ E as well as rela t ivistic corrections and 
focus mainly on therma l SZE as the first trial observation (|Birkinshawlll999l : lEnfilin &: Kaiseiil200Cl : 



Yamada fc Fuiitall2nni 



The study of the SZE has been directed toward understanding the thermal properties of the 
intra cluster m edium (ICM), and it has been detected i n dozens of clu sters of galaxies (see for 
recent reviews, 



Birkinshaw 



19991 : 



Carlstrom et al 



2001 



McKinnon et al 



Rephaelil |2002| ) . Among the numerous 



(199lh first tried to detect the SZE in the 



SZE detections in clusters of galaxies, 
radio lobes with the NRAO 12-m telescope at 90 GHz using the double-subtraction method. They 
observed 4 FR-II radio galaxies whose lobe sizes were smaller than the beam size, and obtained 
upper limits on the antenna temperature fluctuation (~ 0.1 mK, or y ~ 10"^). In this paper, we 
report on the results of the imaging observation of a giant radio galaxy B1358+305, extending much 
further in angular size than the beam size (~ 80" in our observation) using the Nobeyama 45-m 
telescope at 21 GHz. A two dimensional imaging study of a region in a radio galaxy is expected to 
provide the most reliable limit on the SZE in a radio galaxy by resolving the substructure within 
and the emission sources in the field of view. 

The organization of th e paper is as follows . In Section 2 we briefly review the overpressured 
cocoon model (see for detail. lYamada et al ] ll99d ). In Section 3 we describe the features of our target 
B1358-I-305, and the observation and analysis procedures. In Section 4 the results are presented. 
We discuss the possible uncertainties in the estimation of the SZE amplitude in B1358-I-305 and 
their implications for radio galaxy models in section 5. Finally, we summarize our conclusions and 
future prospects in the last section. 
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2. Overview of Overpressured Cocoon Model 



To provide an interpretative fr amework, we review our model of the overpressured cocoon 



(see for details 



Yamada et al 



19991 '). which is an extension of the classical overpressured-cocoon 



model of iBegelman &: Cioffil (jl989l ). In this model, the jets from the AGNs form external shocks 
at the interaction with the surrounding IGM. Due to the action of the reverse shock (or internal 
shock), the jets are compressed. In the FR-II radio galaxies, the compressed jet matter at the 
locations of the internal shock is believed to appear as hotspots. The IGM compressed by the 
external shock is in pressure balance with the internal shock compressed jet matter at a contact 
discontin uity. A high pressure gas clump forms around the AGN, which expands supersonically into 



the I GM (jBegelman fc CiofHll989l : lNathlll995l : iHeinz. Revnolds. Sz Begelmanlll998l : lYamada et al 



19991 ). Hereafter we use the term "cocoon" as the region consisting of both the jet-supplied matter 
and the shock-compressed IGM throughout this paper0 The cocoon expands as its internal energy, 
supplied by the jet, increases with time. In this paper, for simplicity, we first treat is a single free 
parameter that characterizes the cocoon, and assume that the cocoon has an approximately spatially 
uniform thermal electron distribution. Thus the parameter is to be taken as an appropriately 
intensity- weighted mean value over the cocoon, and hereafter we denote the mean value as (.^e)- This 
assumption may not be fully consistent with current observations, but its simplicity can provide a 
basis for feasibility evaluation and interpretation as a first step (see Section [5.21 for the discussion 
of this assumption). 

In these models, most of the matter that has passed through the shock is assumed to be almost 
fully thermalized except for a small fraction, which is assumed to be non-thermal. The density in the 
cocoon is insufficient for effective radiative cooling and, hence, the cocoon is expected to conserve 
nearly all of the kinetic energy of the jets. Furthermore, the cocoon is expected to remain hot for 
a relatively long time due to the low radiative cooling efficiency (an effective life time could be 
tiife ~ 10^ yrs or more: see I Yamada et al .11 19991 ). Even taking into account the PdV w ork against 
the surrounding IGM, the cocoon would remain hot for ~ 10^ yrs ( Yamada et al.lll999l ). Therefore 
the energy contained in a cocoon can be written as 



(3) 



where Ljet is the kinetic luminosity of the jets, and tufe is life time scale of the AGN (typically 
~ lO'^ - 10^ yrs). 

For the thermal SZE as written ([2]) with the thermal pressure in the cocoon and the low value 
of (^e), the total energy of the cocoon becomes C/tot = f^th/(l — iCe)) = / PedV/{l — {^e)) which 
is proportional to y (Eq.[2]). This estimate indicates that SZE, which is proportional to J Pedl, 
provides a good measure of thermalized electrons and the jet-supplied energy contained in the 



Note that the definition of the term "cocoon" can change with different authors, and should not be confused with 
that defined herein. 
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cocoon. We employ this interpretative framework and attempt to measure the energy contained in 
electrons in a radio galaxy. 

Since direct measurement of Lj^t from observation is difficult, equation ([3]) is not used to 
estimate the expected value for the Compton y parameter. Instead, we use the synchrotron emission 
argument and make the assumption that the synchrotron-emitting electrons are generated by diffuse 
shock acceleration with a small acceleration efficiency (again, we assume is the same for the 
internal and external shocks). The minimum energy of the synchrotron-emitting electrons C/syn can 
be estimated from the synchrot ron luminosity Lsyn using the minimum energy condition, which in 
the simplest form (|Moffetlll975l ) is given by 



syn 



1 J 9a / vrc \ 

2 1 T 



1/2 



a+l 
2a + I 



a+l/2 a+1/2- 



4/7 



a+l 



-1 



^3/7^4/7 
-^syn ) 



(4) 



where a is the ratio of energies of non-thermal protons and non-thermal electrons, a is the photon 
index of synchrotron emission, z/„ and z/; are the upper and the lower limits of observation frequen- 
cies, and V is the volume of the synchrotron-emitting region, respectively. Once we obtain the 
minimum energy C/gyn) we can estimate C/tot ^ t4yn/(Ce)j where the inequality takes into account 
the larger internal energy compared with the minimum energy. If we adopt the small acceleration 
efficiency of the shock-accelerated electrons (^e)) the thermal electron energy contained in a cocoon 
is evaluated by Compton y parameter as. 



y 



TJleC- 

rUpC^ V 



Pedl, 



9.2 X 10 



-5 



SOOkpc J \0.05 



V 

-1 



^ rUeC^ V 



U^yJV 



9.7 X 10~i^ergs cm~3 



(5) 
(6) 



along the line of the sight that passes through the center of the cocoon. Numerical values are 
based on the assumption that the axial ratio of the cocoon is 1:3 (Figure [H) and the adopte d 
cosmological parameters {Hq = 72 km s~^ Mpc~^, Qjv/ = 0.3, I^a = 0.7, ISpergel et al.ll2003l ). 
For the mean acceleration efficiency (£p) , we take the small value that is observationally su ggested 



0.05 - 0.1: e.g., ISturner et al.l 119971 : Baring et al 



199S 



Loeb fc WaxmanI 120001 ). The 

synchrotron-emitting electron energy in B1358-I-30 5 was estimated to b e Usyn/V ~ 5.7 — 9.7 x 10~^^ 



ergs cm in the 10 MHz < < 10 GHz band by Parma et al. (jl996l ) (Yamada & Parma, private 
communication), who assumed that the energy ratio of non-thermal protons and electrons, and the 
filling factor of the synchrotron electrons are equal to be unity. Note that equation ^ indicates the 
Comptonization parameter y, or the amplitude of the intensity decrement of SZE is proportional 
to the length within the radio galaxy along the line of sight, so that a giant source can yield a large 
value of y. 
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Observations and Data Reductions 



3.1. Strategy of Search for SZE in a Radio Galaxy 

Our basic plan of observation in search for the SZE signal in a radio galaxy and the feasibility 
evaluation are provided below. Compared to the case of galaxy clusters, the spatial extent of the 
cocoon relative to the rad io lobes is less understood both theoretically and observationally (see 
McNamara &: Nulsenll2007l . and references therein). We, therefore, develop the basic observational 
plan paying a special attention to the extent of the cocoon as follows. 

1) The amplitude of the intensity decrement of SZE is proportional to the length along the line of 
sight (Eq.[6]). Hence a large radio galaxy should be selected as a target. 

2) In order to reduce the cancellation of the intensity decrement of SZE by the emission from 
the radio lobes, the observation frequency should be higher than that used in the observations of 
the synchrotron emission of the radio lobes. Synchrotron emission has a power-law spectral energy 
distribution oc with a typical a about 0.7, and its intensity decreases with v in the GHz regime. 

3) Taking into account the potentially wider spatial extent of the cocoon in comparison to the radio 
lobes that is assumed in the classical hydrodynamic models (e.g.. iBegelman &: Cioffilll989l ). a wider 
region than the radio lobes should be mapped. 

4) Since the expected signal of the SZE is weak and the cocoon structure is not well understood, 
substructures and emission sources within the mapped region should be carefully resolved and 
removed. To resolve the substructures and to take account of the proportionality of the y parameter 
to the length scale along the line of sight (Eq.[6]), a giant radio galaxy whose angular size is also 
large should be selected target. 

I n addition to the point 2, the amplitude of the thermal SZE is known to be maximum at ~ 90 
GHz ()Birkinshawlll999l ). hence the first choice for the observation frequency is 90 GHz. However, the 
beam size decreases with the observation wavelength, and thus the observation frequency cannot be 
too high in order to perform the imaging observation of a wide region (the point 4) in a reasonable 
timescale. The beam size at 90 GHz of the NRO 45-m telescope is about 20", which is too small 
to image the large central region of B 1358+305 (> 5') to achieve a reasonable signal-to- noise ratio 
within a realistic timescale (see discussion on the feasibility in Section 13.1.11 below) . Additionally 
the atmospheric stability becomes worse for higher frequencies in this waveband, or the typical 
system temperature Tgys ~ 250 - 900 K at 80 - 90 GHz is much higher than that at 21 GHz (Tgys w 
100 - 150 K). Hence, the required observational time becomes unrealistically long. 

We chose the single-dish observation as the first trial rather than interferometric observation. 
This follows from the fact th at the hydrodynamic models predict a smooth er distribution of SZE 
compared to galaxy clusters (jBegelman &: CiofRlll989l : lYamada et al.lll999l ). In other words, the 
smooth SZE signal would be resolved out if an interferometric observation were carried out. How- 
ever, interferometric observations would be useful for detecting the SZE for more compact objects 



-7- 



with known structures. In addition, high angular resolution observation with interferometers can 
also be used to remove emission source contaminations within a field of view. Both single-dish and 
interferometric observations have their advantages and disadvantages. In order to avoid the risk of 
resolving out of the possibly extended SZE signal, we use the single-dish observation strategy to 
measure the total flux. 

Given these considerations and the low luminosity at 10 GHz, we decided to observe the central 
region of a giant radio galaxy B1358-I-305 with the NRO 45-m telescope at 21GIIz. The typical 
beam size of the 45-m telescope is ~ 80" at this frequency. The low system temperature (typically 
~ 100 - 150 K) and a wide bandwidth (2 GHz) of the HEMT 22 receiver installed on the 45-m 
telescope led to a decrease in the necessary observation time for an SZE search in a faint giant 
radio galaxy. The detailed set-ups and observation method are described in Section 13.31 below. 



3.1.1. Feasibility Estimation 

The imaging sensitivity for the observation with the 45-m telescope at 21 GHz is about 
2Tsys/rj/y/B = 7.55 mK s^^^, where the typical system temperature Tgys is taken to be 135 K, 
the aperture efficiency r] is estimated to be 0.8, and the bandwidth B is taken to be 2 GHz for one 
of the dual channels, respectively. In order to obtain a 3cr detection of the expected intensity fluc- 
tuation from equation ^ or 0.30 mK with the beam width A^hpbw — 80" in the central 6.7' x 6.7' 
region, the imaging sensitivity estimated above requires an observation time of about 67.3 hrs. 

The estimated observation time implies a deep imaging. As the expected SZE decrement 
in the B1358+305 is small, we examine the cumulat ive contributions f r om ba ckground (or fore- 



ground) discrete sources to the intensity fluctuation. iFranceschini et al.l (|l989l ) calculated the ex- 
pected fluctuations of the antenna temperature (AT^/Ta) as a function of beamwidth at 6 cm, 
produced by randomly distributed sources within a flux range 10~® < S'(Jy) < 10 (see Figure 3 in 



Franceschini et al.lll989l ). We assume that contributing sources have a spectral index a = 0.7 on av- 
erage {ly oc i^~"), and estimate the fluctuation of the brightness temperature ATb/Tb{^ ATa/Ta) 
at 21 GHz to be 1.1 x 10~^. In terms of Compton y parameter, this temperature fluctuation cor- 
responds to y = 4.7 X 10~^ from equation ([1]). The evaluated confusion limit is only marginally 
smaller than the expected SZE (Eq.[6]). Hence in the data analysis procedure, emission sources in 
the field of view should be carefully removed from the integration. 



3.2. The Giant Radio Galaxy B1358-|-305 



B1358+305 is one of the largest rad io galaxies of FR-II type at a moderate redshift {z = 0.206, 
Parma et al.lll996l: ISaripalli et al.lll996l ). Employing cosmological parameters for the A-dominated 
flat universe (jSpergel et al.l 120031 ) . its projected size amounts to ~ 926 kpc, extending over ~ 10' 
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in the plane of the sky ( Parma et al.lll996l )rl It has a relatively low radio luminosity (i-i.4GHz = 
1.9 X lO^^W Hz-i), and the spectral age was estim at ed to about '^2.5 - 7.5 x 10 yrs from the 



multi-band imaging observations (|Parma et al.l 119961 ). iParma et al.l (|l996l ) further calculated the 



advance speed of the shock head of B1358+305 as R/vspec — 0.02c - 0.03c {R is the lobe length, 
Tspec is the spectral age, c is the speed of light, respectively). Prom the balance of the external ram 
pressure and the jet ram pressure, they concluded that the environment density of B 1358+305 is 
quite low (ne < 10^^ cm~^), and B1358+305 is overpressured even for transverse direction, for an 
ambient temperature as high as that of a typical poor cluster environment (~ 10^ K). This makes 
B 1358+305 a good target in the sense that it resembles the classical model. 

B 1358+305 has an FR-II morphology with two aligned lobes with a radio core, which resides 
closer to the northern lobe. As SZE appears as the decrement of the CMB in the frequency range 
ly < 200 GHz, any compensating radio emission should be avoided. Since the lobe intensity rapidly 
decreases towards the region slightly south of the core, we focused on that region in order to avoid 
the lobe emission contamination as much as possible (Figure[l]). 



3.3. Observations 

Observations of B 1358+305 were performed with the HEMT 22 receiver on the Nobeyama 45-m 
telescope from 2001 February 28 to March 18. The receiver has 2 GHz bandwidth, and is equipped 
with dual channels which can simultaneously receive two circular polarization components, which 
enables the reduction of the required observation time. The near-central region in B 1358+305 of 
size 6.7' X 6.7' around {050,650) = (13''58'^25.0^ +30°32'0.0") was raster-scanned with 40" spacing, 
which created a map composed of 11 x 11 pixels. Since the jet axis is nearly along the declination, 
raster-scans were performed along the RA-DEC coordinates in turn. The scan speed is about 25" 
s~^, resulting in the acquisition of one map in ~ 180 s. The exposure time was about 74.8 ksec in 
total, corresponding to 618 s for each pixel. The system noise temperature was about Tgys = 140 K 
on average, ranging from about just below 130 K to nearly 150 K. The noise produced by the 
fluctuations in the atmosphere were reduced by the simultaneous observation of an off-source point 
via the beam switching technique. The off-source point is about 400" away in azimuth direction 
to the east and was simultaneously observed at 15 Hz. The antenna temperature was calibrated 
using the chopper-wheel method, while the optical depths of the atmosphere were measured by 
elevation scans. Since the atmospheric optical depths were about 0.04 throughout the observing 
period, the flux was not corrected for atmospheric absorption. The main beam size ^hpbw was 
estimated to be 81.2" by observing 3C273. With this beam size we can image a large (6.7' x 6.7') 
area in a reasonable timescale, simultaneously resolving the substructures in the mapped region. 
For the pointing and flux calibrator, we used 3C286 about 7.5° away from the center of the image. 



iParma et al.l l|1996l ) estimated the size of B1358+305 to be 1.34 Mpc making use of cosmological parameters 



Hq — 100 km s ^ Mpc ^ and go — 1.0. We update cosmological parameters in this paper. 
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The pointing offset was monitored about every hour, and was typically sm a.ller than &' thr oughout 
the observing period. The flux of 3C286 was assumed to be 2.56 ± 0.2Jy (jOtt et al.lll994l ). and it 
was used to evaluate the aperture efficiency of the telescope. Due to the good weather conditions, 
the aperture efficiency was stable during the latter half of the observation period (rms of the flux 
variation is about 5%). In the former period, the stability was slightly worse (the flux variation 
reached about 10%), and we carefully corrected each observation using the efficiency measured just 
before the observation. 



3.4. Data Reduction and Analysis 

Raster-scans were performed in two orthogonal directions (along the right ascension and the 
declination coordinates) in turn in order to obtain pairs of orthogonally scanned maps taken as 
close together in time as possible. Within a raster-scanned map, antenna temperature fluctuations 
induced by the time variation of the atmospheric conditions, instability of the detectors or other 
instruments and so on, appeared along the scanning direction (scanning noise). The scanning noise 
in the data was reduced by c ombining the two orthogonally scanned ma ps in the Fourier space via 



the basket-weaving method (ISieber. Haslam. &: Salterjll979l : lKundll993l ) 



Although we tried to avoid contamination from lobe emissions into the mapped region, some 
of the scans suffered from lobe emissions in the operation of beam switching. Due to the limited 
throw amplitude (~ 400") and to the rotation of the position angle of the beam switch, emissions 
from one of the lobes entered into the reference held, resulting in over-subtraction in a part of 
the mapped region (less than 1/4 of the mapped area at most). This over-subtraction generated 
artiflcial intensity decrements. We identifled the affected pixels in each scan and re moved them in 



the iii tegration. In the identiflcation process, we made use of the 10.6 GHz map of lSaripalli et al 



(|l996l ) (see Figure[I]), and set the position of the reference beam 400" away from t he main beam. 



Since at 21 GHz both lobes are expected to be dimmer than displayed in the map in lSaripalli et al.l 



(jl996l ) (the spectrum is a decreasing function of frequency) , the procedure above does not overlook 
the affected pixels. Due to this partial masking in the integration, the noise level is not uniform 
in the obtained integrated map, ranging from 0.72 mJy beam^^ to 1.2 mJy beam^^ (Icj level): in 
spite of the loss of some pixels, however, the root mean square of the intensity in the integrated 
map reached 1.22 mJy beam~^ with emission sources. 

In order to evaluate the noise level without emission sources, we constructed two independent 
maps by dividing the maps taken in the former and the latter halves of the observing period, and 
subtracted the latter from the former. In this process, emission from true sources should cancel out 
and the resultant map should contain only random noise. The rms fluctuation of the differential 
map constructed in this way is 1.01 mJy beam~^. 
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Results 



Figure [2] shows the integrated map of the central region of B1358+305. One can easily identify 
the AGN component and the innermost edges of the two radio lobes in this figure, which appear 
also in Figure [TJ On the right edge of the map exists a bright region, with a peak intensity of ~ 6 
mJy beam~^ (source "D"). We failed to find any identified source corresponding to this bright 
spot in any catalogue from the radio to X-ray energy band, though, it indeed appe a rs on the VLA 
NVSS 1.4 GHz map and the Effelsberg 10.6 GHz map obtained by ISaripalh et al.l (|l996l ). Taking 
into account a relatively large beam size (^hpbw = 81.2"), we recognize it as an unknown point 
source. 

Since emissions from the sources smear out the SZ flux decrement, it is necessary to subtract 
them from the map. However, it is quite difficult to obtain the correct zero flux level for mainly 
two reasons. One is related to the larger uncertainties in models of cocoons in comparison with 
clusters of galaxies. In the case of a galaxy cluster, a well-established model obtained by X-ray 
observation (so-called /3 model) is avai lable. Thus, one caii set the baseline (zero flux level) at a 



sufficient distance {9 ^ ^core; see e.g., Komatsu et al.l[200ll ) from the center with the help of the 



X-ray observation. On the other hand, in the case of the theoretical model for the cocoon, there 
are many unc ertainties in its morphology . For example, the density distribution of the surrounding 
medium fe.g.. iKaiser &: Alexandeii 11993 ) or the magnetic fields alter the shape of the cocoon (see 



Burns. Norman. &: Clarke 



1991 



for a review of numerical simulation s, and for examples of MHD 



simulations, see e.g.. lClarke. Norman &: Burnall986l : lLind et al.lll989l l. The other reason is that our 
deep imaging of this region approaches the confusion limit, so the background fluctuation impedes 
the establishment of a deflnitive zero flux level (see Section fS.l.ip . Therefore, we could not ascertain 
the absolute intensity of the sources included in the map at 21 GHz, and it is difficult to correctly 
subtract the emission from the true emission sources. Instead we evaluated the intensity fiuctuation 
in the mapped region, and re-defined the baseline so that the average fiuctuation became zero. In 
order to calculate the intensity fiuctuation, we masked the emission source regions. These regions 
around the emission sources are indicated by white lines in Figure [3l The baseline of the map is 
re-defined after the source removal so as to set the average intensity of the remaining pixels to be 
zero. 



4.1. Structure of the Cocoon 

In order to improve the effective signal-to-noise ratio, we constructed projections of the two- 
dimensional map along two orthogonal directions. We can also expect that these projections refiect 
the structure of the cocoon more clearly than the original two-dimensional map if the data contains 
the SZE signal from the cocoon of B1358-I-305. As can be seen in Figure [U the jet axis is almost 
perpendicular to the right ascension. Thus, we can examine the structure perpendicular (parallel) 
to the jet axis by projecting the map along the declination (right ascension) coordinates. In the 
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process of projection, we average the pixels along each column or row without using the masked 
pixels corresponding to the emission sources that are indicated with the white lines in Figure [3l 

Figured] (a) shows the projection along the jet axis (the declination coordinate). A relatively 
dark region appears at about +2' west of the central column. Comparing Figure H] (a) with the 
original two-dimensional map (FigureE]), one can see the decrement appearing in Figure [H (a) 
mainly picks up the dark region between the unidentified source labeled "D" in the map and the 
AGN (labeled "A"). However, as the statistical errors are large (error bars in Figure H] correspond 
to la) even though they are reduced by ~ 1/ y^Np in averaging over each column {Np is the number 
of pixels included), we cannot rule out the possibility that the projection obtained here is consistent 
with a flat distribution (no SZE signal accompanied 61358+305). As shown in FigurelH B1358+305 
has a relatively large axis ratio, and it is reasonable to approxirnate it as a cylinder with an axis 



along the jet axis. Hydrodynamic simulations (e.g.. iLoken et al.lll992l : iKaiser &: Alexander 



1999 



Scheck et al.ll2002l ) show that the pressure pe is nearly uniform inside a cocoon when the density 
distribution of the surrounding IGM is nearly flat (there will be a difference in at internal and 
external shocks, and accordingly would not be spatially uniform as is suggested by observation, 
but this assumption is discussed separately in Section [5^ below) . Since y (x pel (Eq. [2], where 
I is the length cut through the cocoon along the line of sight), the projection along the jet axis 
is expected to present l{6) in the isobaric cocoon, where 6 is the angular distance projected onto 
the sky from the map center. For a cylindrical cocoon, therefore, the intensity decrement Aly is 
proportional to — \/l — 6'^. We overplot the expected profile oc — Vl — 0"^) based on equation 
([6]) with fiducial values on Figure [4] (a) as a dashed line, offset so that the average intensity is 
equal to zero. The obtained distribution is not inconsistent with the expected one or with the 
flat distribution. This implies the dominance of the statistical errors over the signal reflecting the 
possible cocoon structure in the observed region. We estimate the standard deviation from the flat 
distribution, 

\AIi,\ = 0.28 mJy beam'^ (7) 

which sets the upper limit yupp 

< 1.04 X 10-^ at the 95% confidence level (Eq.[l]). 

We further examine the allowed upper limit for a Compton y parameter by a chi-square test 
of the data for the projection along the jet axis and the expected profile from a cylindrical cocoon 
model AI,y oc — \/l — f^^id that as long as A/j, < 0.32 mJy beam"^ (or y < 1.19 x 10~^), 

these two distributions can agree at the 95% confidence level. Herewith we denote these as |A/,y|'^ 
and yjpp. 

We have also constructed the projection onto the jet axis. If the cocoon has a prolate figure 
and lies in the plane of the sky, the pressure scale length along the jet axis is larger than the 
declination range of the observed field (see Figured]). Hence, the projection onto the jet axis would 
display a more moderate decrement in contrast with the projection along the jet axisH Figure H] 



■^If the jet axis is not straight or the inclination angle is not small, the projected profile will be different from 
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(b) shows the projection onto the jet axis. There clearly exist a peak about 0.7' north of the central 
row and two decrements toward both edges, which are not consistent with the expectation from 
the isobaric cocoon model filled with thermal electrons (Section 2). These decrements are due to 
the two dark regions north and south of the unidentified source "D" (at around a^o ~ 13^58™20'^: 
see Figure[2]). The standard deviation from a flat distribution derived from the projection in this 
direction (onto the jet axis) is =0.44 mJy beam~^, which is larger than the statistical error 

estimated from the noise level of the original two-dimensional map in the process of the projection 
(1.01/y^A'p, ~ 0.27 mJy beam~^ on average). Even though the dark regions appearing in Figured] 
may indeed suggest the existence of SZE signals, we cannot rule out the possibility that these are 
not due to B1358+305. 

In order to examine the dominant source of the fluctuations in FiguresH] (a) and [4] (b) , we 
analyzed the differential map in Section 13.41 in the same manner as the total map. Specifically, the 
differential map was constructed by dividing the observation period into the former half and the 
latter half. FigureO shows the fluctuations in the projections of the differential map overplotted with 
dashed lines on the total map (solid lines, the same as Figured]). The projections on the both axes 
have similar amplitudes (the rms of the projection along the jet axis is 0.27 mJy beam~^ and the 
projection onto the jet axis is 0.43 mJy beam~^) as those of the total map, but their distributions 
do not resemble the total map. Taking into account the fact that real signals, including foreground 
or background emission sources in the field of view, do not appear in the differential map, we 
conclude that the fluctuations in Figurel5]are largely caused by the excess atmospheric noises at 21 
GHz. 

5. Discussion 

Due to the lack of definite information about the position of the edge of the cocoon, it is 
difficult to define a concrete baseline. For this reason we evaluated the relative variation of the 
intensity in the map. In order to improve the signal-to-noise ratio, we made two projections on 
orthogonal directions. A differential map analysis in Section [3.41 shows that the main source of the 
relative fluctuations in these projections is the excess atmospheric noise rather than the true SZE 
signal induced by the hot electrons in B1358+305. 

We discuss possible factors that may reduce the true amplitude of SZE in B1358+305 in 
comparison to the estimate with equation ([6]), including a non- uniform distribution of thermal 
electrons. In addition, the degree to which our results can affect the previous studies on pressure 
balance of B1358+305 (Section 15. 3p is examined. Finally, we direct attention on radio galaxies 
in clusters of galaxies to demonstrate the usefulness of the SZE as a probe of the energetics of 



the one described here. However, iParma et al.l (| 19961 ) concluded that a bend of the jet or the inchnation are not so 
significant, from the contrast of the intensity of northern and southern lobes. 
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radio galaxy in a more general sense to guide future observations. It is to be emphasized that the 
intensity fluctuations discussed in this section should not be confused with a true SZE signal in 
B1358+305, but only as a rough guide for a quantitative discussion. 



5.1. Thermal Evolution of the Cocoon 

The expected amplitude of SZE is estimated using a steady mean acceleration efficiency (^e) 
(Eq.[B]). If (^e) is significantly different from assumed, the expected value of y is also affected. In 
particular, different cooling timescales for the thermal and non-thermal electrons may alter (^g) 
from the fiducial value (=0.05 in Eq.[6]) in time, especially in the central region where the oldest 
population is contained. Since an underestimate of (^e) leads to the overestimate of the expected 
y amplitude (Eq. [6]), we examine whether the different cooling timescales for the thermal and 
non-thermal electrons can significantly alter (^e) using a simple argument about exhaustion of 
the internal energy. Since the cooling timescale is strongly dependent on the electron energy 7, 
we discuss synchrotron electrons of high 7 (7 > 10^ or more) and thermal electrons (7 < 10 at 
Tg < 100 keV) separately. 



Parma et al.l (|l996l ) derived the spectral ages of the non-thermal electrons for both the northern 
and southern lobes using the compiled multi-frequency observation of B1358+305 from 325 MHz 
to 10550 MHz, finding the break frequency to be 2.4 GHz from the integrated spectral energy 
distribution. Further fitting of the multi-frequency images with the energy loss formula including 
synchrotron and CMB inverse Compton losses leads to spectral ages of rgpec for B1358+305 in the 
range 2.5 x 10^ yrs < Tgpec ^ 7.5 x 10^ yrs. Since the synchrotron cooling timescale and inverse 
Compton cooling timescales are proportional to 7^^, lower energy electr o ns ha ve a longer cooling 



timescale. The simple scaling of the above estimation by iParma et al.l (|l996l ) indicates that an 
electron that emits 325 MHz synchrotron light has a spectral age of about a few times of ~ 10^ 
yrs, or even in the extreme case ~ 10^ yrs at best. We adopt 10^ yrs as the typical cooling time 
scale for the non-thermal electrons. 

On the other hand, a relativistic Maxwellian distribution of high temperature thermal electrons, 

7^/3^exp(-7/9) ^ 

Pe{l)dl = ^^(031^^7, (8) 

spans over 7 < 10 for ksT^ ^100 keV (where K2{x) is the second order modified Bessel function, j3 = 
v/c, and G = kBTe/ nieC^ is t he no n-dimensional electron temperature). We estimate the cooling 



timescales following ISarazinI (jl999l ). Inverse Compton scattering of CMB photons, synchrotron 
self-Compton scattering. Coulomb cooling, and bremsstrahlung are taken into accou nt as cooling 
proce sses. The density of the IGM and in B1358+305 is quite low [ue ^ 10~^ cm~^; iParma et al 



19961 ). and thus bremsstrahlung and line cooling are negligibly small compared with the other 
processes. In this low energy regime (7 < 10) and in the very poor environment (the density 
of the surrounding IGM was estimated to be 1.4 x 10^'^ cm~'^ < ?^IGM ^ 8.4 x 10"'^ cm~^; see 
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Table 3 of lParma et al.lll996l ). the major cooling mechanisms are inverse Compton cooling and the 
Coulomb cooling, but the cooling time is longer than the Hubble time (tcqqI ^ 10^^ yrs). Therefore, 
the thermal energy of electrons is conserved for a typical cooling time of the synchrotron electron 
(~ 10^ yrs). Even if we adopt a cooling timescale for electrons that radiate low frequency 325 MHz 
photons (tcooI ^ 10^ yrs), it is unlikely that the cooling time of thermal electrons is lower than 
that of non-thermal electrons. In other words, as long as the radio lobes of B1358+305 are bright 
in synchrotron emission, the thermal energy of electrons in B1358+305 does not change in time. 
Therefore, it is unlikely that the lack of the detection of an intensity decrement in B 1358+305 is 
attributable to the loss of thermal electrons by cooling. 

In addition to the cooling proces ses described above, some models emphasize the importance 
of adiabatic expan sion (llto et al.ll2008l ') or backflows of jet matter that escapes thermalization at the 
shock interaction (^Scheck et al.lbo02|). The se models predict a lower thermal pressure of electrons 
than the model of Begclman & Cioffi (119891 ). Quantitative evaluation of these effects has yet to be 
determined, but it is to be noted that our model tends to overestimate the thermal pressure of the 
electrons in a cocoon. If these effects that reduce the thermal pressure of electrons are significant 
in the evolution of a radio galaxy, much more sensitive observations will be necessary. 



5.2. Non-Thermal Electrons in Radio Lobes 

In this paper, we adopted a model based on the classical overpressured cocoon model with a 
constant and uniform mean shock acceleration efficiency (^e)) and a uniform distribution of thermal 
electrons within a cocoon as the simplest and most optimistic interpretation of the observation 
(Section [2]). However, in spite of the search for X-rays in radio galaxies in galaxy clusters, the clear 
evidence for non-therm al electron accelerati on at the external shock has not been found except for 



objects such as Cen A (jCroston et al.ll2009l ). Although this lack of clear evidence may reflect the 
limited observational sensitivity or angular resolution, the uncertainty in determining (^e) is also 
likely to contribute. Especially, at the internal shock or termination shock of the relativistic jets, 
can be as large as close to unity, while at the non-relativistic external shock is much smaller. 
In this case, the radio lobes within the shocked shell would have only a small population of thermal 
electrons, whereas the shocked shell has a large fraction of thermal electrons as in the first adopted 
model (Section [2]). 

The spatial SZE profile of this model would differ significantly from that adopted (see Section 
14. ip . The amplitude of SZE caused by non-thermal power law electrons depends on the parameters 
describing the electron energy distributions (e.g., the power law index, the lower and upper cut-off 
energies, and the break energy if any). The ratio of thermal and non-thermal SZE amplitudes is 
consequently dependent on these parameters. If the non-thermal SZE amplitude is much smaller 
than thermal SZE in the radio lobes, the SZE spatial profile for this case would not reflect the 
length passing through the cocoon, but the length passing through the shocked IGM shell. There- 
fore, it may have a flatter SZE distribution as compared to that in Section [4. II (oc — Vl — 0"^)- This 
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possibility makes the detection of SZE much more diffi cult (see e.g.. iPfrommer. EnBlin &: Sarazin 
20051 : IColafrancescd 120081 : iHardcastle &: LoonevI l2008l ) . On the contrary, it may also be possible 
that the non-thermal SZE amplitude is similar to the thermal SZE. In this case, the spatial profile 
does not differ from that in Section I4.H and detection of SZE is unlikely to be significantly more 
difficult compared with the cocoon model ffiled with thermal electrons (Section [2]). Multi- frequency 
observation will be necessary to discriminate these distributions (see Section 15.41 below) . Addition- 
ally, it is possi ble that a mature giant radio galaxy does not expand supersonically in the lateral 
direction (e.g., iKonar et al.ll2009l ). In this case, the SZE signal peaks at the head of the external 
shock and decreases toward the central region where we observe (see below). Since the diffuse 
shock acceleration efficiency problem remains to be solved, we cannot exclude either of these sce- 
narios with different acceleration efficiencies. In order to solve this problem and to obtain a deeper 
understanding of the subsequent dynamical and thermal evolution of a radio galaxy, progress in 
both theoretical studies of shock acceleration efficiency and higher sensitivity observations will be 
required. 

Radio images and some X-ray images of radio galaxies suggest that the central region close to 
the nucleus lack the energetic electrons that many hydrodynamic models suggest. If this applies to 
B1358+305, our results would correspond to genuine atmospheric noise. However, the X-ray images 
of a FR-II radio galaxy MS0735+74221 a nd Cygnus A show extended diffuse X-ray emissions that 
fill the c entral regions (see for a revi ew, iMcNamara &: NulsenI |2007| . and references therein) . In 
addition, IHardcastle &: LoonevI (j2008l ) observed bright FR-II galaxies at 90 GHz with the BIMA, 
and found that some of their sample (3C 20 and 3C 388) have extended emissions that cover the 
central regions. Although there are ambiguities in the distributions of the thermal and synchrotron 
electrons, these results imply that some FR-II radio galaxies seem to have electron populations at 
the central regions away from the jets and hotspots. We cannot tell if B1358+305 is deficient of 
electrons because of a small number of past observations or limited sensitivity of our results. Hence 
we conclude that our results are largely due to atmospheric noise, but we cannot completely rule out 
the small SZE signal in B1358+305, and much higher sensitivity is required in future observations. 



5.3. Expansion of B1358-|-305 

Since only upper limits on the y parameter were obtained, we compare our results with past 
studies of B1358+305 in the case where (1) the true pressure is close to that derived from the upper 
limit, and (2) the true pressure is significantly lower than the obtained upper limit. 

Consider the first case. Since the observed region does not extend over the entire cocoon, 
the perpendicular projection onto the jet axis will show a flatter intensity distribution than the 
parallel projection (i.e., the pressure scale length along the jet axis is longer than the declination 
range of the observed field: see Figure[T|). Hence, it is more easily subject to the contamination 
from the background emission sources and the excess atmospheric noise. On the other hand, the 
right ascension range of the observation is expected to be comparable to the cocoon width and the 



- 16 - 



parallel projection along the jet axis is expected to represent the structure of the cocoon (Section 
4.2). Therefore, we employ the upper limit of the y parameter calculated with the projection along 
the jet axis (|A/,y|2o- = 0.56 mJy beam~^, or yupp < 1.04 x 10~^) in the discussion below. If we 
assume that the obtained intensity fluctuation accompanies B1358+305, then the electron pressure 
derives Pe = 2.20 x 10~^^ dyne cm~^ (Eq.[6]) for /=300kpc. If we take values obtained by the 
chi-square test i/upp. the equ i valent electron pressure becomes pe = 2.51 x 10^^^ dyne cm^^. On 
the other hand, iParma et al.l (jl996l ) deduced the jet ram-pressure to be > 1.2 X 10^12 dyne cm , 
which is comparable to the equipartition value under the assumption t h at B1 358+305 expands 
supersonically. This result is quantitatively consistent with iParma et al.l (|l996l ) in the sense that 
B1358+305 expands supersonically against the surrounding IGM. However, the results obtained 
using th e upper limit i n our observation (2.20 - 2.51 x 10^^^ dyne cm ^) are almost double the 

dyne cm~^). This difference would come from the 
and may imply a larger mean acceleration efficiency 



result of Parma et al 



(11993) (> 1.2 X 10 
assumed parameter value (^e) = 0.05 (Eq.[6] 
((^e) ^ 0.1) in B1358+305 and significant contribution from non-thermal SZE of smaller amplitude 
than thermal SZE. 

Since the intensity fluctuations obtained turned out to be dominated by the excess atmospheric 
fluctuations, however, it is also possible that the true thermal electron pressure in B1358+305 is 
much lower than that calculated with y^pp or yupp (the case (2)). Because B1358+305 resides in 
a poor environment, it is difficult to directly measure the IGM pressure around it. Instead we 
estimate the external pressure under the assumption that the environment resembles an extremely 
poor galaxy cluster or a group of galaxies, and assume the IGM temperature to be < 10^ K. For 
the density, we can take the cosmological baryon density at the redshift z = 0.206 as the reference 
value, and describe the IG M density with a density excess parameter e(> 1), njcM = e4.3 x 10~^ 
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cm ^ ( Spergel et al.ll2003l ). which sets the lower limit to the IGM pressure to be > e5.9 x 10 



dyne cm~^. In this case, if the true thermal electron pressure in B1358+305 is less than 10~'^e 
of that inferred from our upper limit, it indicates that thermal electrons do not play a major role 
in supporting the radio lobe against the external pressure. We note that this conclusion strongly 
depends on the estimate of the IGM pressure and the SZE measurement. Deeper X-ray observations 
for the pressure determination of the IGM and much sensitive SZE measurement are required to 
be more conclusive. 



5.4. Comparison with X-ray Observations of Radio Galaxies in Clusters of Galaxies 



Finally, we discuss the X-ray observations in conjunctio n with the SZE in radio galaxies a s 
effective means of probing the pressure of a cocoon (also see iPfrommer. Enijlin &: SarazinI l2005l ) . 
The observations by ROSAT and Chandra satellites have revealed that radio lobes of radio galax- 
ies centered at the cluster of galaxies genera te cavities in the sur f ace brightness distribution of 



the X-ray emitting int r acluster medium (e.g.. iBohringer et al 



McNamara et al 



20001 : 



Fabian et al. 



2000; 



1993; 



McNamara Sz Nulsen 



Carlh. Perlev. Harridll994l : 



20071 ). These are considered to be 
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dueto the pressure in the lob e overwhelming that of the surrounding hot ICM. lHardcastle &: Worrall 
pOOg ) and lLeahv &: Gizanil (2002), however, discovered that the ICM pressure inferred by the X-ray 
emission obtained by ROSAT observation of bright FR-II radio galaxies is much higher than the 
minimum pressure estimated by th e synchrotron luminosit y in the radio band for large (> 100 kpc) 
radio galaxies. A detailed study bv lLeahv fc: Gizanil (j200ll ) confirmed this pressure discrepancy for 
3C388. These results suggest the dominance of the other forms of pressure above that due to the 
synchrotron electrons (and magnetic field) in the radio lobes. 



Hardcastle Sz Worralll (j200d ) argued that plausible candidates for carrying the "invisible" pres- 



sure are non-thermal protons, magnetic field with strength differing from that derived from the 
minimum energy condition (~ B^q), or low energy electrons which do not emit strong radiation in 
the observed energy bands. If we assume that in the internal shock is as small as that in the 
external shock, the radio lobes could have a significant amount of thermal electrons (Section[2]). 
In this case, the pressure of the thermal electrons could be one of the candidates for the origin 
of the "invisible" pressure. Since SZE is sensitive only to electron energy, it is a good probe of 
either the thermal electrons or non-thermal electrons of low energy which do not emit high fre- 
quency emission in the observation bands. The observational different signature between thermal 
and non-thermal SZE is the frequency of zero amplitude (e.g., lBirkinshawlll999l ). Specifically, the 
thermal SZE has the null frequency (= 218 GHz) which is almost independent of the gas tempera- 
ture, whereas the null frequency of non-thermal SZE changes with the electron energy distribution 
parameters. Hence, future multi-frequency observations that covers the null frequency of thermal 
SZE will be a tool to probe the electron contribution to the "invisible" pressure components in the 
cocoon. On the other hand, if of the internal shock is larger than that of the external shock 
and close to unity, resulting in a significant difference in the fractions of thermal and non-thermal 
electrons in the shocked shell and the radio lobes inside, SZE studies become much more difficult. 
In this case, the multi-frequency observation should be carried out over more than one region. For 
example, multi-frequency observations of the head-top region of the bow shock will be able to de- 
termination the contribution of mainly thermal electrons formed at the external shock. This can be 
used to correct the observation of the central region apart from the shock- head, where the thermal 
and non-thermal SZE coexist|l| In summary, the SZE can be a probe of electron contribution to 
the "invisible" pressure in a radio galaxy, but discrimination between thermal and non-thermal 
contributions requires sensitive multi-frequency and multi-region observations in the future. 

For the other sources of pressure, it is unlikely that the magnetic field is significantly different 
from Beq. For example. X-ray observations have shown that majority of the radi o galaxies tend to 
have larger particle energy than that of the magnetic field (for a recent results, see llsobe et al.ll2009l . 
and references therein). Protons, as another possible candidate for the "invisible" pressure, on the 



*For the radio galaxy in a galaxy cluster, another difficult problem is to distinguish the contribution from the 
unshocked ICM. In order to examine this contribution, multi-dimensional modeling of a radio galaxy embedded in 
the ICM is necessary. 
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other hand, would lead to an SZE much smaller than that given by equation ([6l). The existence o f 



high-energy protons should be probed in the hard X-ray and/or 7-ray regime (jScheck et al.ll2002l ) 



6. Conclusion 

We have reported the results of a trial observation of SZE in a giant radio galaxy B1358-I-305 
with the Nobeyama 45-m telescope at 21 GHz. By performing the imaging observation, we have 
obtained the most stringent upper limit achieved for the Compton y parameter in a radio galaxy. 
The obtained upper limit is close to the expected value derived for a low acceleration efficiency 
of synchrotron electrons at the shock (the low value (.^e)) see Eq. [6]), but detailed analysis shows 
that the obtained intensity fluctuation is likely to be caused by the excess atmospheric noise. If 
we assume that the obtained intensity fluctuations in the observed region were due to the thermal 
electrons in the cocoon of B1358-I-305, our results are qualitatively consistent with the supersonic 
expansion of B 1358-1-305 but for the original assumption (^g) being too small. Alternatively, if 
the true pressure of thermal electrons is much lower than that derived from the obtained upper 
limit, our observation is not sensitive enough to derive any definitive conclusions on the pressure 
balance of B1358-I-305. Since the SZE is sensitive only to electron energy, it would serve as a 
probe of the pressure components in radio galaxies as well as galaxy clusters. Future high sensitive 
multi-frequency SZE observations of multi-regions will provide an important information about the 
contributions of thermal or low-energy non-thermal electrons to the total pressure of a radio galaxy. 
This may provide a clue for disentangling the pressure discrepancy between the surrounding ICM 
derived by X-ray observations and the minimum energy of rad io galaxies at the centers of clusters 



of galaxies (jHardcastle &: Worralll 120001 : iLeahv &: Gizanil |2002| ) . Though we have failed to obtain 



a definitive SZE signal with the Nobeyama 45-m telescope, higher frequency observations using 
either a large field-of-view (e.g., > 5'), multi-beam receivers on large (e.g., > 30 m) single dish 
telescopes, or interferometers with a large number of small dishes (e.g., < 50 cm) similar to the 
AMiBA projedH might lead to the detection of SZE of a cocoon. The SZE would provide a new 
observational tool to probe the energetics of radio galaxies, along with the projects to detect low 
freque ncy radio em ission from non-thermal low energy electrons like the Long Wavelength Array 



(LWA: lHarrisll2005l l. 
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^See |http : // amiba . asiaa . sinica . edu . tw/ 1 for the AMiBA project. 
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B1358+305@10.6GHz 
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Fig. 1. — A full image of B1358+305 observed with the Effelsberg 100-m telescope at 10.6 GHz (solid 
contour: reconstructed from an original image in ISaripalli et al.l 119961 ). The region we observed 
(6.7' X 6.7') is indicated by the dashed-line square around the center (indicated by a cross) at 
("50, (^so) = (13'^58™25.0",+30°32'0.0"). A fat cross in the northern lobe indicates the position of 
the radio core. A schematic picture of the cocoon is overlaid as a gray elliptical figure. 
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Fig. 2. — Integrated map of the central region of B1358+305 at 21 GHz. Labels A, B, C and D 
indicate the AGN, north and south lobe edges, and the unknown source (see text), respectively. 
Contours are —3, —2, —1,0, 1, 2, 3 mJy beam~^ (see numbers on the contours). Intensity is normal- 
ized so that the average is zero in the observed region, and the contour levels relative intensities. 
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Fig. 3. — Locations of the emission sources. Masked region corresponding to the AGN, north and 
south lobe edges, and the unknown source on the western side (labeled "D" in the integrated map) 
are indicated with solid white lines. 
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Fig. 4. — Projections along (panel (a)) and onto (panel (b)) the jet axis. Error bar is la. In 
panel (a), relative decrement at about +2' west of the central column is mainly due to the dark 
region north of the source labeled "D" and west of the southern lobe edge (labeled "C"). A dashed 
curve denotes an expected decrement for a cylindrical isobaric cocoon, and a dotted line is the flat 
intensity (no SZE decrement). In panel (b), there is a peak about 0.7' north of the central row and 
two decrements toward both edges. These decrements are due to two dark regions north and south 
of the unidentified source "D" (see Figure [2]) . 
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Fig. 5. — Projections along (panel (a)) and onto (panel (b)) the jet axis of the differential map 
(dashed lines; Section [3.4p . Solid lines are the same as Figured! The differential map projections 
which consist of only random noises have similar amplitudes with the total map, but different 
distribution from it. This indicates the intensity fluctuations are likely to come from the excess 
atmospheric noises (see text). 



